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We study the Abehan flux tube profile of mesonic and the baryonic configurations of 
static quarks below and above the finite-temperature transition in Nf = 2 full QCD. To 
reduce effects of ultra-violet fluctuations we measure Abelian distributions of the action 
density, the color-electric field and the monopole current after fixing to the maximally 
Abelian gauge. Changing the distance R between quarks for fixed T/T^ < 1, one can see 
a clear signal of string breaking for large R. The electric field becomes Coulomb-like and 
the circulating monopole current disappears. We study also the temperature dependence 
of the flux profile for fixed i? ~ 0.7 fm. The disappearance of the squeezed flux is observed 
clearly above T^. Similar behaviors are observed both in the mesonic and the baryonic 
systems. 



1. Introduction 

The study of various distributions inside the flux tubes formed in meson and baryon 
systems is important to understand the internal structure of mesons and baryons and to 
learn more about the QCD vacuum properties. In full QCD the flux tube, which is usu- 
ally modeled by the bosonic string, is expected to disappear (string breaking) when the 
distance between quarks becomes large enough. How the flux tube evolves and how it dis- 
appears when the distance between quarks or the temperature increase beyond respective 
critical values are still open questions to be addressed in this paper. 

In zero temperature full QCD no clear evidence for string breaking has been found 
so far [in El IS]- This is due to poor overlap of the Wilson loop with the two-meson 
channels [Hj. At zero temperature the flux tube has been extensively studied in quenched 
QCD [ Ej while there were only few lattice studies [ El El EI in full QCD. In the latter 
papers the calculations were done in the Abelian projected theory after fixing to the 
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Table 1 

Parameters and statistics of the simulation, together with the integrated autocorrelation 
time Tint- The length of the trajectory is r = 0.25. 



maximally Abelian gauge (MAG). The use of Abelian projected quantities allowed us to 
reduce ultra-violet fluctuations and improve the signal to noise ratio in comparison with 
the non- Abelian quantities. In [O] the profile of the Abelian flux tube between a quark 
and an anti-quark of the length i? ~ 1.6 fm was presented. The three quark system 
has been studied in [ Ej . It has been shown that the baryonic flux tube has a shape 
consistent with the so-called Y-type flux tube in agreement with other numerical results 
[ IH] , and theoretical predictions [ 13 HDl UH 112] • In finite temperature full QCD string 
breaking is much easier to observe using Polyakov loop correlators [ • The state 

created by Polyakov loops has considerable overlap with both the string and the broken 
string states. 

In this paper, we study flux tube profiles in the mesonic and the baryonic systems in 
finite temperature lattice QCD with dynamical fermions at distances below and above 
the string breaking distance. In [ we have studied string breaking in this theory using 
Abelian Polyakov loop correlators. In that paper the transition temperature was also 
determined. Following our previous work [ [71 [Ej we investigate the structure of the 
flux tube with the help of Abelian observables after fixing to MAG. In particular, we 
study the monopole and the photon parts of Abelian flux tube profiles. This permits 
us to get good signal for both mesonic and baryonic flux tubes of large size and to find 
further confirmation of the dual superconductor mechanism of confinement. 



2. Simulation details 



To study QCD with dynamical quarks we consider Nf = 2 flavors of degenerate 
quarks, using the Wilson gauge field action and non-perturbatively 0{a) improved Wil- 
son fermions [ ISj. Configurations are generated on 16^ x 8 lattice at (3 = 5.2, 0.1330 < 
K < 0.1360, corresponding to temperatures below and above the transition temperature 
Te = 213(10) MeV {nt = 0.1344) and at ^ = 5.25 for k < [114|. Details of the sim- 
ulations can be found in [ |2] • Table ^ shows the number of configurations used in this 
study. We perform the MAG fixing on generated configurations employing the simulated 
annealing algorithm [El- The Abelian projection procedure [ E] defines the diagonal 
link matrices 



diag{ul{s),ul{s),ul{s)}, 



(1) 

(2) 
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where ^JJJ(s) can in turn be decomposed into monopole (singular) and photon (regular) 
parts [[IHllinil2ni: 

o;{s) = e';^°^'-{s) + ef'-{s). (3) 

Monopole and photon parts are defined up to irrelevant terms by following expressions: 
^~(.) = -27rY.Dis-s')dlnt^{s) , (4) 

s' 

Of'^'is) = -T.D{s-s')d'X,is), (5) 

s' 

where 

= d,ei{s)~d,ei{s)-2'nn';^{s), e;,is) e n^M = 0,±1,±2, (6) 

du' and dy are backward and forward lattice derivatives, respectively, and D{s) denotes 
the lattice Coulomb propagator. 

We study the Abelian action density 

Pab(.) = f EEC0S(^"^.(^))' (7) 

the Abelian color-electric field 

E^{s)=i9';,{s), (8) 
and the monopole current 

klis) = -^t,y,^dJlAs + fi) . (9) 

We consider three types of Polyakov loops to create static sources: 
Abelian: 

Ll^{s) = exp ^ i E ei{s, t) , L,^{s) = - E Ll^{s) ; (10) 



t=i 



monopole: 



^mon(s) = exp iAY.^T'''\s,t) \ , Lmon(s) = -E^mon(^ 



photon: 



f Aft ^ 1 3 

LU^ = exp |iEC'^(^"'>^)| ' V(^1 = 3 E^ph(^1 • 



(12) 



The vacuum averages of our observables are defined for the mesonic case by 
, . (Pab(s)pQQ(r)) 

{pMs))qq = ,^ (Pab(s)), (13) 



4 




Figure 1. Illustration of a global rotational transformation used for off-axis distances in 
the QQ system. 
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and for the baryonic case by 

(pab(s))3Q = ,^ ,^ (Pab(s)), 



(14) 
(15) 



(E,(S)) 



3Q 



(fcj(s))3Q 

where 



(^3Q(rY)) 

(^3Q(rY)) 

(jr E |£abc|fc,"(g)L°(3i)L''(s2)i^'^(g3)) 
(^3Q(rY)) 



(16) 
(17) 



a=l 



^3Q(rY) = ^ E kafecl L^(5'i)L^(5'2)^'=(5'3) 



(19) 

(20) 



Si denotes the position of the r quark, r is the distance between quark and anti-quark 



r = |si - S2I 



(21) 



and ty is the minimal Y-type distance between the three quarks, i.e. the sum of the 
distances from the three quarks to the Fermat point. 



4 = ^E4 + 2V3^A. 



(22) 
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Figure 2. Expectation values of the non-Abelian ("NA" x25), the Abehan ("Ab" x 3.3), 
the monopole ("Mo") and the photon ("Ph") Polyakov loops (left) and the static potential 
from monopole Polyakov loop correlators (right) for l3 = 5.2. 



Tij = \si — Sj\, S/\ is the area of the corresponding triangle. Eq. (j^ defines ry when all 
angles in the three quarks triangle are less than 27r/3. If one of the angles is equal or 
larger than 27r/3, then [El '"y = J2i>j^ij ~ maxrjj. L°'{s) stands for one of three types 
of Polyakov loops introduced above. The lattice distance is converted into the physical 
scale using the lattice spacing a(/9), i.e., R = ra and Ry = rya. 

For mesonic system we need both on-axis and off-axis distances. In the latter case the 
X-axis was rotated along the line connecting the sources as shown in Figure H 

3. Mesonic system 

3.1. The finite temperature transition and the static potential 

Some features of the finite-temperature full QCD have been studied in our previous 
paper [Ej. We show expectation values of various Polyakov loops in Figure |21 (left). 
One can see that all Polyakov loops are smooth functions of the temperature and that 
Abelian and monopole Polyakov loop behaviors are qualitatively the same as that of the 
non-Abelian Polyakov loop, while the photon Polyakov loop is almost constant across the 
transition [E]- The transition is a smooth one (cross-over) in contrast with the first-order 
phase transition in quenched QCD. 

The heavy quark potential V{r, T) is determined, up to an entropy contribution, from 
the Polyakov loop correlator: 

^V{r,T) = -\n{L{s)L\s')) . (23) 

In [ E| we measured non-Abelian, Abelian, monopole and photon potentials using for 
L(s) respective Polyakov loop operators. When r oo, 

{L{s)L\s')) ^ \{L)\\ (24) 
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Figure 3. The profile of tlie action density of tlie mesonic system at T/T^ = 0.94. The 
sources are made of Abehan (left), monopole (center) and photon (right) Polyakov loops. 
The QQ distances are R/ro = 0.91 (top), 1.59 (middle), 2.25 (bottom). 



where \{L)\'^ ^ 0, since the global Z3 is broken by fermions. 

In Figure 121 (right) the monopole part of the heavy quark potential is shown. We have 
analyzed the static potential assuming that at temperatures T < the Polyakov loop 
correlator can be described in terms of the string and two heavy-light meson states. Then 

{L{s)L\s')) = e-(^o(^)+^-i-(-'^))/^ + e-2^(^)/^ , (25) 

where for the non-Abelian potential V^tring('", 7") is given in Ref. [ E], while for the 
monopole potential we used the linear term only. The energy E{T) can be written as 

E{T) = ^Vo + m{T), (26) 

where m{T) can be considered as a constituent quark mass [122], Vq is a self-energy of the 
static source. The second term of Eq ()26p gives rise to the flattening of the potential at 
large distances and describes the effect of creation of a pair of constituent quarks. 
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Figure 4. The profile of tlie color-electric field for the mesonic flux tube at T/T^ = 0.94. 
The sources and distances are as in Figure El 



3.2. Profiles of the mesonic system 

As we discussed in the Introduction, in full QCD the flux tube disappears at large dis- 
tances due to the light quark- ant i- quark pair creation. In this subsection we demonstrate 
how this happens for the mesonic flux tube when the Q — Q distance R is increased at 
fixed temperature T. We measured profiles of the action density, the color-electric field 
and the monopole current at T/Tc = 0.94 for R/vq = 0.91, 1.59 and 2.25. The static 
quark- ant i- quark pair was created by Abelian, monopole, or photon Polyakov loops. 

Figure 121 shows profiles of the action density. We show both the two-dimensional contour 
plot and the three-dimensional plot for every distance. In the contour plots, the horizontal 
axis is always along the line connecting the static sources. First two columns of this figure 
indicate that the flux tube persists for distances R/tq = 0.91, 1.59 while for R/tq = 2.25 
it seems to disappear, leaving only lumps of the action density around the sources. Note 
that the latter distance is roughly the distance where the static potential, shown in Figure 
121 becomes flat. The difference between two columns is that in the monopole part the 
data has much lower noise and two-peaks at locations of the quark and anti-quark are 
much less pronounced. In the photon part only two lumps around sources can be seen at 
all distances. 

In Figure m we show the profile of the color-electric field in the plane including the QQ 
axis, whereas the azimuthal component of the monopole current is plotted in Figure 
The figures suggest that the color-electric field is squeezed into a flux tube due to the 
magnetic supercurrent circulating around the QQ axis and creating the solenoidal electric 
field. The agreement with this picture has been found before in SU(2) gluodynamics and 
in both quenched and full QCD at zero temperature [ |2S1 1211 IHl HI- Figures El and El 
present first evidence that this picture is also correct in full QCD at finite temperature [ 
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Figure 5. The profile of tlie monopole current for tlie mesonic flux tube at T/T^ 
The sources and distances are as in Figure El 



0.94. 



IT^ . One can see from Figure lU that for two smaller distances the monopole part of the 
electric field has the form of a sourceless solenoidal field [0]. This field changes direction 
at some distance away from the QQ axis and cancels in this region the Coulomb field 
created by the photon part. As a result the total Abelian electric field has the form of 
a flux tube. Comparison of the monopole part of the electric field for R/ro = 0.91 and 
R/tq = 1.59 reveals that the electric field becomes weaker with increasing distance. It 
virtually disappears for R/tq = 2.25, thus providing clear evidence for the string breaking 
at R/tq = 2.25. String breaking at this distance can also be seen from the monopole part 
of the magnetic current, depicted in Figure |S1 Indeed the azimuthal monopole current 
disappears at this distance, while it is clearly seen at smaller distances. It is known from 
studies of the SU(2) theory that the magnetic current calculated from monopole sources 
reproduces almost exactly that calculated with Abelian sources. This can also be seen 
from the first row of Figure El In the second row the Abelian part shows only random 
fluctuations which is due to stronger noise for this case in comparison with the monopole 
part. For the photon part, as expected, we see only noise at all distances. 
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Figure 6. The profile of tlie monopole part of the action density (top), the color-electric 
field (middle) and the monopole current (bottom) at (left to right) T/T^ = 0.80, 0.94, 
1.00 and 1.28 for R/ro ~ 1.36. 



3.3. Temperature dependence of the mesonic system 

It is also interesting to see how the flux tube of given length changes with the temper- 
ature. We measured the same observable as in the previous subsection varying now the 
temperature, while keeping the QQ distance R approximately constant. 

We show in Figure IHl the profile of the monopole part only because the signal for it is 
much better than that for the Abelian part. For T/Tc = 0.80 and 0.94 we see the clear flux 
tube in all our observables. We can see also that all observables decrease with increasing 
temperature. This behavior is similar to the flux tube evolution with increasing distance 
discussed in the previous subsection. 

At the critical temperature T = Tc, there seems to remain a touch of a flux tube as seen 
from the action density, although signals of the electric field and the monopole current 
are very weak. Note that the existence of a weak flux tube at this temperature would be 
in agreement with the continuous character of the transition, which is expected to be a 
crossover. 

In the deconfinement phase at T/T^ = 1.28, no tube profile is observed. This is consis- 
tent with the disappearance of the linear potential at T > T^. 

4. Baryonic system 

4.1. Static potential 

Now we discuss the system constructed from three static quarks. We measure the static 
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Figure 7. The static three quark potential at T/Tc = 0.94. 
three-quark potential defined as: 

V^Q{st,n,st,T) = -^ln(^Y.l^}^abc\L^(^i)L\n)L%n)) ■ (27) 

Since the baryonic fiux tube has Y-shape at T = [ 13 HI], we use ry defined by 
Eq. (jUJ as a typical distance describing the three-quark system. 

Note that the Y-shape of the fiux tube is also supported by results for the static 
potential in quenched QCD [E]^- As is seen from Figure [7| the static three quark potential 
has a tendency to fiatten at large ry. At finite temperature not only the lowest energy 
state (baryon) but also states of broken string, shown in Fig. |H1 may contribute to the 
potential (ITfl) . 

The breaking of the string occurs due to creation of pairs made from a light quark and 
a light anti-quark, qq, from the vacuum. Contrary to the mesonic case discussed in the 
previous section, breaking of the baryonic string may happen at seven different stages, 
depicted schematically in Figure |H1 

1. Process QQQ + qq QQq + Qq. 

The most energetically favorable string breaking corresponds to creation of one 
virtual meson, qq, from the vacuum. The light quark from this pair and two heavy 
quarks form a baryon system, QQq. The chromoelectric string inside the QQq- 
baryon is spanned between the two heavy quarks. The light anti-quark and the 
third heavy quark form a heavy-light, Qq, meson state. There are three possible 
ways to break the Y-string by one light meson, they are shown in the second row of 
Figure |H1 

2. Process QQQ + 2qq Qqq + 2Qq. 

Another possibility is the creation of two light meson states, 2qq, from the vacuum. 

^However the data of ref. [|2S] give preference to A type potential 
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Figure 8. Terms contributing to the correlation function of three Polyakov loops ()27jl . 
The large circles represent the static quarks (Q), the small circles correspond to the 
dynamical quarks (g). 



Then the string becomes completely broken and the breaking produces two heavy- 
light mesons, Qq and one baryon made from the heavy quark and the two light 
quarks, Qqq. There are also three possible combinations of the final states depicted 
in the third row of Figure |H1 

3. Process QQQ + 3qq — > qqq + 3Qq. 

Finally, three light mesons can be created from the vacuum producing three heavy- 
light mesons Qq and one ggg-baryon made from the light quarks (the bottom of 
Figure |HI) • 

All these states must be taken into account in the correlator of the three static quarks. 
Thus instead of the two-exponential formula ()25|) we will have the five-exponential ansatz 
which includes also the original unbroken string state (the first term below): 



_|_ ^-^Qq/T-mQQq/T 
_|_ 3 g-2(mQq+mQqg)/T 
_|_ g-3(mQ9+m,,,)/r| ^ 



(2J 



This simplified formula takes into account only the linear potentials between the heavy 
quarks as well as the broken string states, disregarding other possible interactions between 
the quarks. It is known [ ESI HEl 0] that both in quenched and unquenched QCD the 
monopole part of the mesonic potential is linear at large as well as at small distances. In 
the study of the mesonic potential in full QCD at finite temperature [[Hj we found again 
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Figure 9. The monopole contribution to tensions of chromoelectric strings inside the 
QQQ and QQq states for /? = 5.2 and /? = 5.25. We show for comparison the tension 
of the string obtained in the meson state [QQ] [ 14^, and the shaded area shows the 
quenched value of the QQ string tension [[29j. 

that the monopole part of the potential can be fitted well by Eq. ()25|) which includes only 
the linear potential in Vstring{r,T). Thus, we may expect that the correlator (^Tj) of the 
monopole Polyakov loops (jllj) can be described by the simple form of Eq. (j28|) . 

In Eq. (j28p we explicitly used two different string tensions, one for the tension of the 
Y-shaped unbroken string state (the first term in the r.h.s. of ((211)) and another for the 
string tension of the QQq baryon. On the classical level there is no difference between 
strings inside QQQ and QQq states. The only difference is in the shape of strings : the 
QQQ state has Y-shape string while in the QQq state the chromoelectric string between 
heavy quarks forming an approximated straight line. In the QQq state the sea quark, 
q, should be located somewhere inside the string. However, on the quantum level the 
sea quark may interact with the string and the string tension may, in general, get an 
additional correction from the quantum motion of the sea quark. 

In Eq.()28|) the dependence on the self-energy Vq is factorized similarly to the mesonic 
case (I25|) . The expression in the curly brackets contains terms depending on the con- 
stituent masses of light quarks in the Qq, QQq, Qqq and qqq states {mQQg, mQq, mqqq 
and ruqqq, respectively). The potential depends not only on Ry, but also on Rij 
and this is the reason of the broad distribution of points in Fig. [7| Note that the two 
completely broken string states, Qqq + 2Qq and qqq + ?>Qq (corresponding to the last 
two rows in Figure |H1) provide distance-independent contributions to the potential (j2Hl)- 
This independence exists only for the monopole part of the potential, which neglects the 
Coulomb interactions between quarks, and, consequently, van der Vaals forces between 
the mesons and the baryons. Thus, for the non-Abelian potential or for the corresponding 
Abehan part of it the contribution of the Qqq + 2Qq and qqq + 3Qq states to the potential 
should be dependent on the distances between the quarks. 

We have analyzed the data for the three quark correlation function following Ref. [ 
IS], but with the ansatz ()28p which takes into account the breaking of the string. The 
results for both string tensions are presented in Figure M We show for comparison the 
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Figure 10. The same as in Figure IHl but for the average of the effective quark masses in 
the Qq meson and the QQq baryon, (mg^ + mQQq)/2. We also show the effective quark 
masse of the Qq meson obtained from two-point Polyakov loop correlators in Ref. [E]- 



tension of the string obtained in the mesonic case in Nf = 2 QCD Ref. [ E] and in 
the corresponding quenched value [ One can see that string tensions are almost 
independent of the temperature in the region close to the phase transition. Moreover, the 
string tension of the QQQ system is slightly higher than the string tension of the QQq 
state. Thus, we conclude that the sea quark lowers the energy of the string. 

Another interesting result presented in Figure M is that the tension of the unbroken 
string in the baryon, agqq, is slightly lower than the tension of the string in meson, (Tqq. 

Using the fitting function ()28p we obtain the sum of the constituent masses of the light 
quarks in the Qq and QQq states. As one can see from the form of the fitting function ()28|) 
we can not get the values of these masses (as well as the other mass parameters) sepa- 
rately using the monopole contribution to the correlation function of the three quarks. 
In Figure ITUl we present the averaged constituent quark mass, {rriQg + mQQg)/2, which 
shows a tendency to drop with increasing temperature. We also depict in this Figure the 
constituent quark mass in the Qq meson which was obtained from two-point Polyakov 
loop correlators in Ref. [|14J. Our data shows that the constituent quark mass in the 
QQq-haijon is higher than in the mass in the heavy-light meson, Qq. 

4.2. Profiles of the baryonic system 

In Fig^^we show profiles of the action density and the color-electric field in the baryonic 
system. For the static quark sources we use the monopole Polyakov loops. For Ry/vq = 
1.96, 2.61 the linear part of the static potential exists and the profile of the action density 
has Y-shape. The electric field is also squeezed into Y-shape. However, the noise is much 
bigger than that for the T = case [ El E] • 

We are using the correlator of three Polyakov loops to describe the baryonic system at 
finite temperature. Thus there is no a Y-shape junction in the operator itself, contrary 
to the case of the Wilson loop used at zero temperature [El E|- Nevertheless we clearly 
see that the flux distribution has a Y-shape geometry with a junction, similar to the zero 
temperature case [El EI- 
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Figure 11. The action density (top) and the color-electric field (bottom) of the baryonic 
system. The sources are made of the monopole part of Polyakov loops. The distances 
are Ry/rQ = 1.96(left), 2.61(center) and 3.26(right), and the temperature is fixed to 
T/Tc = 0.80 («; = 0.1330). 

For large enough Ry string breaking is expected. We see this effect in the action 
density in Fig^Jat Ry/fQ = 3.26. The color-electric field is too noisy though to observe 
string breaking. In Fig^J we show the color-electric field and the monopole current 
distributions for T/T^ = 0.80 and Ry/r^ = 1.96. The electric field is squeezed into a 
Y-shape fiux tube as we already discussed. One can see circulating monopole currents 
around the color-electric field in each slice. In the plane where the color-electric field is 
subdivided into two parts, the circulating monopole current is not a simple circle anymore. 
This indicates the possibility for the formation of two circulating currents if the distance 
between quarks is large enough. 

In Fig^Jwe plot the profile of the baryonic system at various temperatures for Ry/r^ ~ 
2.0, monopole Polyakov loops represent static quarks. The behavior of the fiux tube profile 
is similar to the meson case. When the temperature increases, the density of the action 
and the electric field decrease and almost disappear in the deconfinement region. 

5. Conclusions 

We have studied the Abelian profiles of the action density, the electric field and the 
monopole current in the mesonic and the baryonic systems below and above the finite 
temperature transition in Nj = 2 full QCD on A^j = 8 lattices. For the source we used 
Abelian, monopole and photon Polyakov loop operators. In the mesonic case, we see a 
clear flux tube profile for some, not too large, quark- ant i- quark distances R. However, 
even in the confinement phase, fiux tube breaking is observed for large enough R, in 
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Figure 12. The profile of tfie color-electric field (left) and the monopole current (right) 
at T/Tc = 0.87 in a three-quark system. 




Figure 13. The profile of the monopole part in the baryonic system at (left to right) 
T/Tc = 0.80, 0.94, 1.00 and 1.28. These figures correspond to physical distance Ry/ro ~ 
2.0. The top row shows the action density, the bottom row shows the color-electric field. 
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agreement with expectations. We have also studied the temperature dependence of the 
mesonic flux tube profile keeping the distance R fixed at about 0.7 fm. We found that 
the value of the flux density decreases smoothly as the temperature approaches T^. In the 
deconfinement phase, the flux tube disappears. But at the critical temperature, Tc, the 
flux tube still exists. 

In the baryonic case, we have calculated the static three-quark potential in finite tem- 
perature full QCD. However, our statistics is not good enough to study the baryonic 
string breaking in detail. Nevertheless the flattening of the baryonic potential is actually 
observed for large Ry- The flux tube profile and the string breaking phenomenon are 
consistent with the meson case. It is interesting that we could observe the Y-shape profile 
without introducing any junction in the baryon creation operator. 

We are performing now the calculations on 24^ x 10 lattice to check for finite volume 
effects. 

Acknowledgements 

This work is supported by the SR8000 Supercomputer Project of High Energy Accelera- 
tor Research Organization (KEK). A part of numerical measurements has been done using 
NEC SX-5 at RCNP of Osaka University. T.S. is partially supported by JSPS Grant-in- 
Aid for Scientific Research on Priority Areas No. 13135210 and (B) No. 15340073. The 
Moscow group is partially supported by RFBR grants 02-02-17308, 01-02-17456, DFG- 
RFBR RUS 113 1739/0, RFBR-DFG 03-02-04016, grants INTAS-00-00111 DFG-RFBR 
436 RUS 113/739/0, and CRDF awards RPI-2364-MO-02 and MO-011-0. 

REFERENCES 

1. S. Aoki et al [CP-PACS Collaboration], Nucl. Phys. Proc. Suppl. 73 (1999) 216. 

2. B. Bolder et al, Phys. Rev. D 63 (2001) 074504 

3. DIK Collaboration (V.G. Bornyakov et al.), |hep-lat /0310011 ; V. Bornyakov et al., 
Nucl. Phys. Proc. Suppl. 106 (2002) 634. 

4. C. W. Bernard et al, Phys. Rev. D 64 (2001) 074509. 

5. G. S. Bali, K. Schilling and C. Schhchter, Phys. Rev. D 51 (1995) 5165; R. W. Hay- 
maker, V. Singh, Y. C. Peng and J. Wosiek, Phys. Rev. D 53 (1996) 389. 

6. H. Ichie, V. Bornyakov, T. Streuer and G. Schierholz, Nucl. Phys. Proc. Suppl. 119 
(2003) 751; Nucl. Phys. A 721 (2003) 899. 

7. V. G. Bornyakov et al [DIK Collaboration], arXiv:hep-lat/0401026 . 

8. T. T. Takahashi, H. Suganuma, Y. Nemoto, H. Matsufuru, Phys. Rev. D65 (2002) 
114509. 

9. S. Kamizawa, Y. Matsubara, H. Shiba and T. Suzuki, Nucl. Phys. B 389 (1993) 563. 

10. M. N. Chernodub and D. A. Komarov, JETP Lett. 68 (1998) 117. 

11. D. S. Kuzmenko and Y. A. Simonov, Phys. Lett. B 494 (2000) 81. 

12. Y. Koma, E.-M. Ilgenfritz, T. Suzuki and H. Toki Phys. Rev. D64 (2001) 014015 

13. C. DeTar, O. Kaczmarek, F. Karsch and E. Laermann, Phys. Rev. D 59 (1999) 031501. 

14. V. G. Bornyakov et al [DIK Collaboration], arXiv:hep-lat /0401014j , V. Bornyakov 
et al, Nucl. Phys. Proc. Suppl. 119 (2003) 703. 

15. S. Booth et al [QCDSF-UKQCD collaboration], Phys. Lett. B 519 (2001) 229. 



17 



16. G. S. Bali, V. G. Bornyakov, M. MuUer-Preussker and K. Schilling, Phys. Rev. D 
54 (1996) 2863 

17. A. S. Kronfeld, M. L. Laursen, G. Schierholz and U. J. Wiese, Phys. Lett. B198, 
516 (1987). 

18. T. Suzuki et al, Phys. Lett. B347 (1995) 375; [Erratum-ibid. B 351 (1995) 603]. 

19. J. Smit and A. van der Sijs, Nucl. Phys. B355 (1991) 603. 

20. T. Suzuki, S. Ilyar, Y. Matsubara, T. Okude and K. Yotsuji, Phys. Lett. B347 
(1995) 375; Erratum: ibid. B351 (1995) 603. 

21. M. Gao, Phys. Rev. D 40 (1989) 2708; see also: Ph. de Forcrand, G. Schierholz, H. 
Schneid er and M. Teper, Ph ys. Lett. B 160 (1985) 137. 

22. H. Satz, arXiv:hep-ph/0111265 

23. Y. Koma, M. Koma, E.-M. Ilgenfritz, T. Suzuki, M.I. Polikarpov, Phys. Rev. D68 
(2003) 094018, 

24. Y. Koma, M. Koma, E.-M. Ilgenfritz, T. Suzuki, Phys. Rev. D68 (2003) 114504. 

25. C. Alexandrou, Ph. de Forcrand, A. Tsapalis, Phys. Rev. D66 (2002) 094503. 

26. H. Shiba and T. Suzuki, Phys. Lett. B 333 (1994) 461; J. D. Stack, S. D. Neiman 
and R. J. Wensley, Phys. Rev. D 50 (1994) 3399. 

27. T. Suzuki and I. Yotsuyanagi, Phys. Rev. D 42 (1990) 4257. 

28. M. Llischer, P. Weisz, JHEP 0207 (2002) 049. 

29. O. Kaczmarek, F. Karsch, E. Laermann and M. Lutgemeier, Phys. Rev. D 62 (2000) 
034021. 



